Autocorrelations and Thermal Fragility of Anyonic Loops 
in Topologically Quantum Ordered Systems 
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Are systems that display Topological Quantum Order (TQO), and have a gap to excitations, hard- 
ware fault-tolerant at finite temperatures? We show that in models that display low d- dimensional 
Gauge-Like Symmetries, such as Kitaev's and its generalizations, the expectation value of topo- 
logical symmetry operators vanishes at any non-zero temperature, a phenomenon that we coined 
thermal fragility. The autocorrelation time for the non-local topological quantities in these sys- 
tems may remain finite even in the thermodynamic limit. We provide explicit expressions for the 
autocorrelation functions in Kitaev's Toric code model. If temperatures far below the gap may 
be achieved then these autocorrelation times, albeit finite, can be made large. The physical en- 
gine behind the loss of correlations at large spatial and/or temporal distance is the proliferation of 
topological defects at any finite temperature as a result of a dimensional reduction. This raises an 
important question: How may we best quantify the degree of protection of quantum information in 
a topologically ordered system at finite temperature? 

PACS numbers: 05.30.-d, 03.67.Pp, 05.30.Pr, 11.15.-q 



I. INTRODUCTION 

The perseverance of information over long times in the 
simplest of memory devices is related to the existence of 
large associated autocorrelation times. The storage of 
information is intimately tied to the breaking of ergodic- 
ity at scales much smaller than the autocorrelation time. 
Classical information can be reliably stored in magnet- 
ically or in electrically (permanently) polarized materi- 
als. From the physicist's perspective, this reliability is 
directly linked to the existence of an order parameter (its 
macroscopic magnetization or polarization) which char- 
acterizes a collective and robust property of the material 
below an ordering transition temperature. At its core, 
non-ergodicity implies the existence of a generalized order 
parameter (e.g. the overlap parameter of spin glasses). 

The reliable storage of quantum information is a real 
challenge. The uncontrolled interactions between a quan- 
tum system and its environment or measurement appara- 
tus introduce noise (errors) in the system leading to de- 
coherence of pure quantum superposition states. Fortu- 
nately, quantum states can, in principle, be encoded fault 
tolerantly and be protected against decoherence, thus pre- 
venting loss of information [1| . This idea lies at the heart 
of TQO systems as first advanced by Kitaev Q • Assum- 
ing that errors are of a local nature, topological quantum 
memories (e.g. surface codes Q) seem to be intrinsi- 
cally stable because of physical fault-tolerance to weak 
quasi-local perturbations. However, are these quantum 
memories robust to thermal effects? 

In this work, we analyze the effect of temperature on 
zero-temperature (T = 0) topologically ordered quantum 
systems 0,01' sucn as Kitaev's Toric Code [|| and Honey- 
comb models Q and generalizations thereof. To this end 
we need to present two concepts that were introduced for 
the first time in our previous work [|| . One is the concept 



of finite-T Topological Quantum Order (TQO), and the 
other of rank-n TQO. In that same work we studied the 
thermal fragility of topological operators in D = 2 lattice 
models. Our results [6| concerning the singular character 
of the T = TQO in one notable system (Kitaev's Toric 
code model) have later been reaffirmed in work by Castel- 
novo and Chamon Q in their study of the topological en- 
tanglement entropy. In the present work we will present 
extensions of our ideas to higher spatial dimensions D 
and expand on the physical reasons leading to thermal 
fragility. In particular, we show that a general gauge 
theory in D spatial dimensions in a system with periodic 
boundary conditions displays rank-n = k D TQO. Never- 
theless, although a thermodynamic phase transition may 
occur, the system is thermally fragile. We investigate not 
only the thermodynamic but also the dynamical aspects 
of thermal fragility, and in cases such as Kitaev's Toric 
Code model wc also obtain exact analytic time-dependent 
results thanks to our duality mappings [f|. 



II. LANDAU ORDERS VS TQO 

Before defining TQO, and to put this latter concept 
in perspective, let us briefly review the rudiments of a 
Landau order parameter. The Landau order parameter 
is customarily associated with the breaking of a global 
symmetry. The existence of an order parameter, a macro- 
scopic property measuring the degree of order in a state of 
matter, is directly associated to the phenomenon of spon- 
taneous symmetry breaking (SSB). This concept, that 
involves an infinite number of degrees of freedom, is so 
fundamental to condensed matter and particle physics 
that many excellent textbooks (see, for example, [8[) have 
spent entire chapters (or even a full book |9() describing 
it. For the present purposes, wc illustrate the concept in 
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the simple case of a ferromagnct. A piece of iron at high 
temperatures it is in a disordered paramagnetic phase. 
Below a certain temperature T c the system orders, i.e. it 
magnetizes, and with the appearance of the order param- 
eter (magnetization) there is a breaking of the rotational 
symmetry [Ioj |. In the (ferro)magnetic phase there is a 
net magnetization M that persists all the way to zero 
temperature (where it attains its maximal value). The 
magnetization can, mathematically, be written as a linear 
combination of quasi-local operators V (e.g. V is the local 
spin operator in a Heisenberg model). The main point 
to stress here is that the operator V may distinguish be- 
tween different ground states (GSs) \g a ) and \gp) of the 
material 

(9a\V\g a ) ± (gp\V\gp), (1) 

and equivalently at finite temperatures, 

(V) a * (V)p. (2) 

Here and throughout the angular brackets refer to a 
thermodynamic average: for any quantiry- say A we 
have that (A) — tr [pA]. The density matrix p = 
Z- r ey^[-H/{k B T)] with Z = tr [exp[-/3i?]] where /3 
is the inverse temperature. Similarly, we define a (trace- 
class) density matrix p a = Z~ x exp[— H a / [ksT)] (with 
Z a = tr [exp[— f3H a \] the partition function) correspond- 
ing to the Hamiltonian H endowed with terms which fa- 
vor order in the state \g a )- These subscripts (a) are those 
appearing in the thermodynamic averages of Eq.([2]). A 
particular realization of p a for the problems that will in- 
terest us will be given below (Eq. pTj) ). 

Symmetry plays a key role in dictating the fundamen- 
tal properties of matter. Symmetry often generally im- 
plies the existence of conserved charges with unique phys- 
ical consequences. Most Landau orders are inherently 
tied to broken global symmetries. There are symme- 
tries other than global. For instance, the symmetries 
in gauge theories are local; such local symmetries cannot 
be broken [ll|. Recently, a general classification (and 
their physical consequences) of these and other types of 
symmetries was proposed |6|. A symmetry is termed a 
(i-dimensional Gauge-Like-Symmetry (d-GLS) if the min- 
imal non-empty set of symmetry operations operate on a 
d-dimensional spatial volume [y, [l2j. Thus, global sym- 
metries - those of the usual Landau-type - correspond to 
d = D (here the symmetry operators act non-trivially 
on the entire Z?-dimcnsional system), and gauge symme- 
tries (which are local in nature) correspond to d = (as 
the symmetry operators act non-trivially only on quasi- 
local (or d = O-dimensional) regions). General symme- 
tries may lie anywhere in between these two extremes: 
< d < D. The groups associated with such symmetries 
can be denoted as Gd [1, 03 ■ (I n the following we will 
only consider unitary representations.) The statements 
that we will make below pertain to general systems, both 
in the continuum and on lattices. For explicit forms, in 
what follows, we will often provide expressions and refer 



to systems defined on D-dimcnsional hypcrcubic lattices 
of size L x L x • • • x L. 

Let us now define TQO. Given a set of n orthonormal 
GSs {\g a )}a=i,— ,n, w hh 1 < n < N g where N g is the 
total number of GSs of a given Hamiltonian H , T = 
rank-n TQO exists iff for any bounded operator V with 
compact support (i.e. any quasi-local operator V), 

(9a\V\gp) = V S a p + c a p, (3) 

where v is a constant and c a p is a correction that vanishes 
in the thermodynamic limit. This is indeed a condition 
on non-distinguishability of GSs through local measure- 
ments. Here and throughout, we will employ Greek let- 
ters a and (3 to denote orthogonal states in the GS man- 
ifold. Note that Eq. ([3]) applies only to systems with 
degenerate GSs. Following standard conventions, (3 will 
also be employed for the inverse temperature f/(/ssT) 
where ks is the Boltzmann constant. 

General error detection in TQO systems, motivated by 
quantum error detection conditions elsewhere, is given by 

[PoVPo , f M ] = P [V, %]P = 0, (4) 

where Pq = J2 a \9a)(9a\ is the protected subspace, and 
T^'s represent the logical operators which are not part of 
the code 's stabilizer [13| . These operators are non-trivial 
symmetries of H, i.e. 

[H,f fl ]=0, (5) 

and encode the braiding operations that ensure topo- 
logical degeneracy of the GS manifold. In the anyonic 
schemes, these operators represent braiding opeations. 
We will reserve the use of the Greek indices p and v to 
the operators {T^}. It is important to emphasize that the 
T = TQO quantum error detection condition, Eq. 
applies for systems with degenerate and non-degenerate 
GSs (unlike Eq. ([3]) which only applies to systems with 
degenerate GSs). 

Clearly, when c a p = condition ^ implies (HJ). To 
see this, we write the commutator of Eq. (0| longhand 
to find that 

Y,(g a \V\gp)l\g a )(gp\,f,} = v[P ,f^} = (6) 

a, f3 

identically (regardless of the specific symmetry operator 
T M ). This follows from Eq. (0) whenever Eq. ^ holds 
with c a=j Lp = 0. Thus, condition ^ is sufficient (but not 
necessary) to ensure the general error detection condition 
ofEq. 

A finite-T (T > 0) generalization of TQO is provided 
by the condition Q 

(V) a = tr [ Pa V] = v + c aa {L), Va (7) 

with c aa a correction that tends to zero in the thermo- 
dynamic limit. TQO systems @ satisfy both the T = 
as well as the finite-T conditions of Eqs. ^ and (J7J). 



3 



Motivated by the conditions for quantum error detec- 
tion one can propose an extension of (UJ) to finite temper- 
atures. The finite temperature error detection condition 
that we will focus on in this work [l4| is 

[^V^.T^^^rj/^o, (8) 

for all quasi-local operators V. Here, the propagation of 
a local error at finite temperatures (V) causes no harm 
to the logical operators {T M }. In the "typical finite tem- 
perature subspace", the local errors (V = p 1 / 2 Vp 1 / 2 ) do 
not alter the algebra of the operators {T^}. What we will 
ultimately measure at thermal equilibrium are objects of 
the form tr [pViT^V^T,,...]. [For finite times (in which 
equilibration has not set in yet), we will measure finite 
time correlations of a similar form.] If the appearance of 
the local operators Vi does not alter the algebra of the 
symmetry operators T^. Eq.© is the simplest caricature 
ensuring such invariances. In anyonic schemes, the alge- 
bra of such non-local operators- the algebra of {T^}- is 
what enables quantum memories. Similarly the algebra 
of related non-local operators encodes the braiding oper- 
ators that may perform topological quantum computing. 
In the Appendix, we will show that Eq. ([8]) generally can- 
not be satisfied for any system of a finite size nor, more 
generally, in any other system which does not display a 
finite-T transition. The existence of a finite-T phase tran- 
sition is a necessary but not sufficient condition for Eq. 
© to hold. We will now briefly relate a weaker version of 
the finite-T detection condition of Eq. ([5]) to our earlier 
finite-T TQO condition of Eq. 0. This will suggest that 
although finite-T transitions (singularities in the free en- 
ergy) are mandated to ensure finite temperature error 
detection, these transitions cannot be accompanied by 
SSB. Thus, these transitions may be more akin to those 
in gauge theories. To this end, we note that if Eq. ([8]) 
holds then, in particular, the finite-T expectation value 

(V) = (T^yft). (9) 

We now consider an extension of Eq. © that is valid in 
the thermodynamic limit for all quasi-local V 

(V) a = (f^Vfl) a , (10) 

stating that there is no SSB of the symmetries spanned by 
{T M } - at least insofar as any local observable V can de- 
tect. It is worth emphasizing that in Eq. (fT0"|) . the indices 
a and p generally need not be the same. The absence of 
SSB detectable by local observables is the physical con- 
tent of Eq. with 

p a = Z^eM-PiH + KT a )], (11) 

where h a — > 0. 

In our recent work Q, we further developed a sym- 
metry principle for TQO. We related certain symmetry 
transformations of a system to the existence of TQO as 
defined by Eqs. ([3j [7]). We emphasized the fundamen- 
tal role ci-GLSs play in establishing that order and dis- 
allowing SSB of local quantities. We basically proved 



sufficient symmetry conditions for a system to be topo- 
logically quantum ordered: When in a gapped system of 
finite interaction range and strength, the GSs (each of 
which can be chosen by the application of an infinitesimal 
field) may be linked by discrete d < 1 or by continuous 
d < 2 GLSs U € Qd, then a system that satisfies the 
T = conditions of Eq. (0) exhibits finite-T TQO [in 
the sense of Eqs. 00)/. We refer the reader to Q for 
a comprehensive explanation. The quantum error detec- 
tion conditions of Eq. ([8]) are far more restrictive than 
the finite-T TQO conditions regarding the robustness of 
the system to all quasi-local perturbations V . 

Although it is possible in many cases [f| to satisfy the 
weaker version of Eq. (|8|) [that is Eq. (fT0|) ] by the use of 
e?-GLSs, unless they are biased by hand, the logical oper- 
ators {T M } always have a vanishing expectation value at 
any long-time equilibrium finite-T state (i.e. any possible 
Gibbs state) in a system with finite range interactions 

(fja = 0. (12) 

A related, more practical, consequence is that {T M } 
may generally exhibit finite autocorrelation times. That 
is, at all positive temperatures, 

%(*) = <f M (0)f^)) (13) 

with Gf (t) — > as (|t|/r) — > oo; there is an inher- 
ently finite autocorrelation time r, whose size is limited 
by thermal fluctuations (but not by system size). This 
autocorrelation time remains finite even in the thermody- 
namic limit. An indefinitely self-correcting TQO surface 
code can only exist at exactly zero temperature. How- 
ever, it may be that by setting parameters we can tune 
r to be very large. 

In Section ITTTAl we will show how this is explicitly real- 
ized in the Kitacv's Toric code model [Eqs. (j3"9"lH0IH3")) ]. 
Before embarking on an analysis of specific cases, let us 
first analyze general relations. 

In several specific cases, such as Kitaev's Toric code 
model, which we will analyze below, it is possible to find a 
new representation in which an initial local Hamiltonian 
remains local yet the non-local topological anyonic loops 
become objects of low effective dimensionality. For in- 
stance, in Kitaev's Toric code model, the non-local Toric 
cycle loops become point (d = 0) fields. In such cases 
of low effective dimensionality, the absence of finite tem- 
perature SSB in low dimensions guarantees that Eq. (p~2|) 
holds even without performing more details expansions 
or bounds. 

Below, we prove this result for the particular case of 
a symmetry operator which is independent not only of 
the code's stabilizer but also of the Hamiltonian itself. 
Later on, we will show how this follows also when the 
symmetry operator (or, in fact, any non-local operator) 
is not independent of the the arguments on which the 
Hamiltonian depends. When the symmetry operator is 
independent of the Hamiltonian, a transformation exists 
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which turns the symmetry into a d — O-dimcnsional 
operator. Similar to Eq. (|ll[) . we may define 



Z(h, 



tr 



exp[-P(H - h^)} 



(14) 



with no summation over repeated indices (//) implicit. 
If the logical operator is independent of the argument 
of a local Hamiltonian H = J2i Hi (here {Hi} are local 
operators) , then the partition function of Eq. (fT4"|) simply 
factorizes 



ZQi^) = Z x zQi/t). 
From Eq. ([T5j) . the expectation value 



(15) 



Id _., , 
hm ; — - — — Z(n a ) 



= lim 



1 







z(K) = o. 



(16) 



The expectation value evaluated with z(h^) that encom- 
passes only one site, cannot exhibit SSB (in formal terms, 
z(h^) is the partition function of a d = 0-dimensional sys- 
tem). Thus, quite universally, the expectation value of 
any such logical operator vanishes. We emphasize that 
this holds for all systems (both in the thermodynamic 
limit and finite size systems). All that matters in the 
derivation of Eqs. (fT5l [TO)) is that T M is independent 
of the variables {Hi} that are added to form the code's 
Hamiltonian. In Section[TTT]we will work out these and re- 
lated expectation values in detail for Kitaev's Toric code 
model. 

Equation (|12[) may be extended more generally to non- 
local high-dimensional operators (R a ) that need not (i) 
lie outside the code's stabilizer or (ii) be symmetries of 
H . For instance, in lattice gauge theories, R a can pertain 
to a Wilson loop of a divergent perimeter. [l5[ By per- 
forming a low temperature series expansion about the or- 
dered state, one generally finds that, similar to Eq. (fT2")l . 
all non-local operators {R a } of dimension d > have a 
vanishing expectation value at finite temperatures, 



(Ra) = 0. 



(17) 



Equation (|17p follows from an asymptotic perimeter law 
type bound (see, e.g. [lj|) 



\(R a )\<Ae- 



(18) 



with m = 0(L d ) the number of local fields that lie in the 
support of R a and A, c positive constants. [TtJ As seen 
from Eq. (fT5|) . for m — ► oo (as befits any non-local R a ), 
this expectation value vanishes. Hand in hand, 



{Ra(0)R a (t)} ltHoo ^0. 



(19) 



as 



Related results (for both the commutator in Eq. ^ M 
well as the commutator between topological quantities 



are afforded by simple extensions (carried in the Ap- 
pendix) of the Lieb-Robinson bounds known to apply for 
local quantities in spin systems with local interactions 

Consider next d = symmetries, i.e. operators {T M } 
that span only a zero-dimensional volume (or a finite 
number of points on a lattice). For the operators {T M } 
to realize a non-trivial ray representation which leads to 
a topological degeneracy, these local symmetry opera- 
tors cannot commute with one another. Elitzur's the- 
orem [Til ] states that any quantity which is not invari- 
ant under all local symmetries must vanish at any finite 
temperature. Thus, for any quasi-local T M symmetry op- 
erator (including all operators which may be defined on 
any finite size lattice), there can never be a SSB of T M 
and once again Eq. (|12p follows. Considerations similar 
to those of Eq. lfTo]) can be enacted. [T|| Couched in the 
language more commonly used by researchers in anyonic 
quantum computing, we can say that in this case- the 
case of general local (d = 0) operators (which includes 
any finite lattice as a special realization of d = sym- 
metries which are here enforced by the limited physical 
extent D = of the system)- if we given a set of sym- 
metry group operators {T M }, i.e. [H, T M ] = 0, such that 
they form a non-Abclian group or an Abelian group with 
a ray (non- vector) representation, then we can prove that 
(T M ) = 0. For example, in the case of a ray representation 
TpT v = e^">f v T^ (f„) = tr [f v pf}f„] = (f^) = 0. 
From the dimension of the irreducible representation we 
infer the degeneracy of the GS subspace. (We remind the 
reader that for a continuous connected group of symme- 
tries all finite-dimensional ray representations arc equiv- 
alent to vector representations.) 

We conclude this section with a general remark and 
reiterate one of our earlier comments regarding autocor- 
relation times: Although indefinite quantum error detec- 
tion may be ruled out by a system that violates the sim- 
plest finite T error detection condition of Eq.©) or for 
which Eq. (|12p holds, our results do not rule out quantum 
error correction over time scales that can be made quite 
large (albeit still finite) by a judicial choice of parameters. 
In the next section, we will show how although Eq. (fT2")) is 
realized in thermal equilibrium, the autocorrelation time 
may be made large at very low temperatures. 



III. THERMAL FRAGILITY: WORKED-OUT 
EXAMPLES 

In earlier work the presence of a gap in the en- 
ergy spectrum together with the existence of TQO as 
defined above (Eq. ^ with c Q/ 3 = 0) were suggested to 
be sufficient to guarantee the protection of quantum in- 
formation. The physical intuition behind this was that 
properties of the protected (GS) subspace are stable with 
respect to weak local perturbations. How does tempera- 
ture affect this conclusion? Is there any other additional 
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FIG. 1: Elementary (cell) plaquette B p and star A a interac- 
tion terms in Kitaev's Toric Code model. Hollow circles in 
the bonds (links) represent an S = 1/2 degree of freedom, 
while thick (dashed or solid) lines represent the topological 
symmetry operators (see text). 



requirement needed for protection? We just mentioned 
that the simple generalization of the quantum error de- 
tection condition. Eq. ((8|), may generally fail at all tem- 
peratures T =/= 0. It is commonly believed that the exis- 
tence of a finite gap between the ground and first excited 
states protects properties associated with T = TQO up 
to a finite energy scale UbT smaller than the gap A since 
thermal fluctuations are suppressed by the Boltzmann 
factor cxp[— A/ksT]. Our results [|| showed that this as- 
sumption is. in general, incorrect for long times. We orig- 
inally coined the term thermal fragility to describe this 
state of affairs @. We showed that in some prominent 
TQO models, the entropic weight associated with defects 
outweights their Boltzmann penalty of exp[— A/fe^T] at 
any finite temperature: the equilibrium states are always 
disordered. In what follows, we expand on these con- 
cepts, determine the equilibration time in one solvable 
case, and consider extensions of previous examples to 
higher spatial dimensions D. In Section IIII A[ we de- 
rive the autocorrelation function for Kitaev's Toric code 
model. We then review [in Section IIII B| a new high- 
dimensional extension of Kitaev's Toric code model and 
prove that this system displays TQO and a finite auto- 
correlation time. Finally, in Section lill CI we discuss the 
situation for Kitaev's Honeycomb Model. 



A. D = 2 Kitaev's Toric Code Model 

For the sake of clarity and because of historical impor- 
tance we start by summarizing the main findings of Ref. 
regarding Kitaev's Toric Code model Q in D = 2. 



The model is defined on a square lattice with L X L = N s 
sites, where on each bond (or link) (ij) it is defined a 
5=1/2 degree of freedom indicated by a Pauli matrix 
, thus defining a Hilbcrt space of dimension 2 2Ns (see 
Fig. [I]). The model Hamiltonian is given by 

H k = -Y,A s -J2 b p (2°) 



with Hcrmitian operators (whose eigenvalues are ±1) 



A., 



n 



(ij')estar(s) 



n <v 

('ij)Gplaquette(p) 



(21) 



and afj (k = x, y, z) representing Pauli matrices. B p and 
A s describe the plaquette (or face) and star (or vertex) 
operators, respectively, with (Vs, s',p,p') 



[A,, A*] = [B p ,B p ,] = [A s ,B p ]=0, 



(22) 



thus generating an Abclian group called code 's stabilizer 
p|. In the presence of periodic boundary conditions, the 
plaquette and star operators satisfy the constraint 



and the two d 

Zl,2 = 



l[A s =l[B p = l, 

S p 

1 Z2 symmetries are given by 



(23) 



n 



X 



1.2 



n 



(y")eCi, 2 



(24) 



where Ci(C' 2 ) are horizontal and C2(C{) vertical closed 
contours (i.e. loops on the lattice(dual lattice)). The 
logical operators and X\ t 2 commute with the code's 
stabilizer but are not part of it, thus acting non-trivially 
on the two encoded Toric code qubits. 

As shown in Ref. Q Hk is related to Wen's plaque- 
tte model [|[ and to two Ising chains by exact duality 
mappings. Therefore, these three models share the same 
spectrum. The GS (protected subspace of the code) is 
4-fold degenerate (Abelian Z2 X Z2 symmetry) and there 
is a gap to excitations. The spectrum is basically that 
of two uncoupled circular Ising chains (2N S is the total 
number of links of the original D = 2 lattice) 



H K - - / ] gp^p+i - 



E< 

8 = 1 



'S + 1! 



(25) 



with GS energy Eq = —2N S and a gap to the first ex- 
cited state equals to 4 (this value of 4 and not 2 arises as 
in the presence of periodic boundary conditions, only an 
even number of domain walls are possible). This exam- 
ple clearly illustrates the fact that TQO is a property of 
states and not of the Hamiltonian spectrum [6j . 

The elementary excitations of H K are of two types @ 



i**(r)>= n 

(y)er 

i*x(n>= n 

(ii)er' 



l*o> 



:5«(T)|* >, 

s*(r')|* >, 



(26) 
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where r(r') is an open string on the lattice (dual lat- 
tice) and |^o) is a GS. (If T(T') would be closed contours 
which circumscribe an entire Toric cycle then the string 
operators S x ' z would become the Toric symmetries of Eq. 
([24]) .) In the case of the open contours of Eq. ([26]) . the 
operators S x ' z generate excitations at the end points of 
these strings (thus always coming in pairs) with Abclian 
Fractional Statistics (anyons). Excitations living on the 
vertices represent electric charges while the ones living 
on the plaquettes are magnetic vortices. These magnetic 
and electric type excitations obey fusion rules that enable 
Abelian quantum computation. Due to the exact equiv- 
alence between Kitaev's model and the Ising chains, no 
non-trivial finite temperature SSB or other transitions 
can take place. The spectrum exhibits a multitude of 
low-energy states. At any finite temperature, no matter 
how small, entropic contributions to the free energy over- 
whelm energy penalties and lead to a free energy which 
is everywhere analytic jf|. 

The operators of Eqs. (f2"o| are not symmetries of H 
and, for divergent loop size, are specific examples of the 
non-local operators {R a } that we considered earlier (see, 
e.g. Eq. (fT5|) ). In the one-dimensional Ising duality map- 
ping of Eq. (|25|) j w e may represent the string operators 
of Eq. as the creation operators for domain walls in 
the D = 1 Ising model. When T and V intersect at any 
even number of bonds, a representation is 



S Z (T) 



n 



s x {r') 



n 



(27) 



si <s<s 



Pi <P<P2 



In Eq. (|27j) . Si and S2 denote the endpoints of the string 
r of Eq. (J2H). Similarly, p\ i2 are the plaquettes which 
form the endpoints of the string V on the dual lattice. 
When r and V share an odd number of bonds, we can 
represent the string operators as follows 



S Z (T) 



P2 



n 



s x {r') 



n 



(28) 



Sl <S<S2 



Pl<P<P2 



For any finite lattice and T > 



(Z a ) = 



tr [exp[-f3H K ]Z,j,] _ tr [A M exp[- 

tr [cxp[-0H K ]] 
tr [expl-fiHrfXpZuXn] 



tr [exp[-PH K ]] 



tr [exp[-0H K ]] 
= -{ZJ = 0, (29) 



where we have used the property [Hk,Xu\ = = 
\Hk , Zp] and the cyclic invariance of the trace tr , which 
is performed over the eigenstates of afj with eigenvalues 
±1. Similarly, (X^) = 0. Indeed, this is a special case of 
the more general argument depicted in Section ITT1 

This does not, in principle, preclude SSB in the ther- 
modynamic limit. One needs to restrict the configura- 
tions over which the trace is performed. In the thermo- 
dynamic limit, derivatives of the partition function and 
associated free energy need not be analytic single- valued 
functions (when SSB occurs the expectation values de- 
pend on how the limit h — > is taken). To this end, let 



us define the (generating) partition function of the model 
with the constraint (f2"3"|) 

Z = tr [exp[-/3(H K - ^ {h x ^X^ + h z ,M)]] (30) 

ju=l,2 

= [(2cosh/3) JV « + (2sinh/3) Ars ] 2 cosh/3/ii cosh/3/i 2 , 



where h h 



h X4l + h Z4i . Thus, the free energy per 



bond, F = —jj^-\nZ[h K ^ = 0], is analytic for all finite 
/?, and displays a singularity at T = (inherited from the 
Ising chain). That means that no finite-T phase transi- 
tion occurs in Kitaev's model. Moreover, from Eq. (|31|) . 
we can compute the expectation values of the topological 
operators with the result 



(Z. 



lim 



d 



InZ 



Z . j I 



ii*,^0+ d{fih x ii ) " ii«.l-*0+ h 



(Xu) = lim 



d 



InZ 



lim 

h 



lim 

c . u ^o+ h,. 



(Zt) = (Z 2 ) = (Xi) = (X 2 ) = 0. 



t&nh(j3hn), 



tanh( / 3ft, Al ), 



(31) 



This indicates that the existence of a gap in this system 
may not protect a finite expectation value of the Toric 
code operators X\^ or Z\^ for any finite temperature 
T > |, [2(|. At T = these expectation values are 
finite and equal to unity, reflecting the non-analyticity of 
F at T = 0. The physical reason behind this result is 
the proliferation of topological defects (solitons) at any 
finite T. The Boltzmann suppression becomes ineffective 
at sufficiently long times and this might be bad news 
for a robust quantum memory pil ]. In the presence of 
additional fields (sec last terms in Eq. (|3"Tj) ). with the 
two component vector = {h x ^,h z ^) at site /i, we 
may define the two component vector 



1 



(32) 



and set ([H K - £)u=i 2 (h x ,^X li + h z ^Z^), Q M ] = 0) 



Qi = Xin x i + Ztfizi, 

Q 2 = X 2 n x2 + Z 2 n z2 



(33) 



to be the counterparts of two single Ising spins a x c and 
a z d which are located at site numbers c and d of the two 
respective Ising chains (that of the s and that of the p 
varieties) . c and d can be chosen to be any integers such 
that 1 < c, d < N s . These two spins (along with the 
spins appearing in Eq. (|25j) satisfy precisely the same 
algebra and set of constraints as the original variables 
in Kitaev's model [Eqs. <£M HQ 122)] ■ The vanishing 
expectation values of Eq. (|3~lj) for both the finite and 
infinite system can be understood as the statement that 
(°~sc) = (°~pd) = f° r an y one-dimensional Ising system 
when the single on-site magnetic fields h c d — * 0. 

The entropy associated with d = 1 type domain walls 
is logarithmic in the system size N s . By contrast, the 
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energy penalty for these domain walls is finite and size 
independent. As a result, for sufficiently large systems, 
cntropic gains will outweigh energy penalties. In partic- 
ular, in the thermodynamic limit, there is no SSB of the 
d = 1 GLSs at any temperature T > 0. In we showed 
how all static correlation functions may be computed via 
our mapping to the Ising chain. 

We can now address the dynamical aspects of ther- 
mal fragility in Kitaev's Toric code model. From its 
mapping to two uncoupled Ising chains, we can imme- 
diately determine the time autocorrelation functions of 
the Toric code operators. As we have shown, Kitaev's 
Toric code operators {Xip, Zip} map into single spins 
in an Ising chain. Thus, we can employ the results ob- 
tained in |22l . [23[ concerning autocorrelation of spins in 
Ising chains in a system with Glauber-type dynamics. In 
the Appendices (Eqs. (|D3[|D4[) ). we will outline standard 
master equations which used to determine the dynamics. 
When these equations depend only on the system's spec- 
tra (and not the precise real-space form of the GSs) then 
we can use our mapping of Eq. (|33|) to relate the dy- 
namics of the non-local topological quantities of Kitaev's 
Toric code model to single spin dynamics in Ising chains. 
We find that 



obeys, for t > 0, the relation 

dG x , 



(34) 



at 



tanh" (3-1 

X 9 h 

tanh 2 (3 + 1 



t tanh /3 tanh 2/3 



.(35) 



Here, x ls a constant setting the time scale for the evolu- 
tion of the system and Iq is the modified Bessel function. 
For low temperatures, at short times, [23| 



1*1 «X"\ 
the autocorrelation is given by 

Q v . tt . (f) ~ e -v/2[l-tanh2/3] X |t| 

At intermediate times, [23[ 



X 1 < |*| < X 1 



1 



1 - tanh 2/3' 



(36) 



(37) 



(38) 



the autocorrelation of the topological string operators is 
well approximated by a Cole-Davidson form [24j | in the 
frequency domain and stretched exponential in time, 



(*Xn ;intermediate-time (*) 

with the equilibration time 



4x(l- tanh 2/3) 



(39) 



(40) 



The constant r is independent of the system size (for 
large systems) and is finite for all temperatures T > 0. 



If the temperature can be made much smaller than the 
gap, i.e. if (3 3> 1, then r, albeit being finite, can be made 
large (t ~ (77/(8%) exp[4/3]). We recall that in Kitaev's 
Toric code model on the torus (the system with periodic 
boundary conditions) the gap between the GS and lowest 
excited energy levels is equal to four: A = 4 [see the 
discussion after Eq. (|25[) ]. (Similarly, the gap for the 
model on an open surface (open boundary conditions) is 
given by A = 2.) Thus, for Kitaev's Toric code model 
on a torus, in the limit of small temperatures the inverse 
Botzmann factor exp[/3A] scales in the same fashion as 
the equilibration time r. The important feature of the 
scaling of the autocorrelation time in the Kitaev's model 
is that, for large systems, it is system size independent 
at all temperatures. Similarly, in the long-time limit, 



|t|»x" 1 - 



1 



1 - tanh 2/3 ' 

the autocorrelation function is well approximated by 



(41) 



Cjf M ;long-time(*) 

/ 1 + tanh 2(3 
2nx\t\(l- tanh 2/3) 



(42) 



-X(l-tanh2^)|t 



Identical relations hold for the autocorrelators Gz tL (t) = 
(Z^(Q)Z^(t)). Thus, as we increase the system size the 
error rate will always be finite at any T > 0. In other 
words, any topological quantum memory at finite (non- 
zero) temperatures might not sustain self-correction in- 
definitely. 

We now discuss Eq. ([8]). It is readily seen that this 
condition is violated here. For instance, if we choose V = 
afj, and T M = Z\ with (if) £ C\ we have, for any finite 
/3, that 



,% 1/2 /o 



(43) 



as {o~fj,o~ij} = 0. This is, of course, a particular realiza- 
tion of the general result of Appendix [XJ The Kitaev 's 
Toric code model has a T — transition and thus cannot 
satisfy the equilibrium thermal detection condition of Ap- 
pendix\^\ It cannot be ruled out however that Kitaev's 
Toric code model may nevertheless satisfy quantum error 
detection at finite temperature up to a small discrepancy 
nor that error detection may work well for finite time 
intervals. 

On the other hand, by choosing V = of), which has 
the property that 

(1 + B p )af j (l + B p ) = 0, (44) 

and realizing that the GS projection operator is given by 



(45) 



the T = detection condition, Eq. is trivially satis- 
fied. Similarly, if one chooses V = erf) since 



(1 + A s )o-?J1 + A s )=0, 



(46) 



and, in general, for Kitaev's model Eq. (|4|) is satisfied. 
That the T = conditions of Eq. (0| are satisfied also 
follows from the fact that Eq. (J3j implies that Eq. |(3]), 
and as we proved in Ref. [|| Kitaev's model, satisfies the 
T = conditions of Eq. ©. 



Generalizations of Kitaev's Toric Code Model 



There is obviously no connection between the lack of 
thermodynamic phase transition (as in Kitaev's model) 
and the existence of TQO as defined in the Introduc- 
tion. Indeed, we will now explore two models that have 
TQO yet display finite-T phase transitions as signaled by 
non-analyticities in the free energy F. The first model 
we will consider is a D = 3 extension of Kitaev's model 
and the second a D = 3 gauge theory. In both of 
these systems, no SSB occurs (and TQO is indeed mate- 
rialized). Nevertheless, the system's free energy displays 
singularities at finite temperatures. 

The D = 3 Kitaev's extended model (KE3Z?) in a cubic 
lattice has a Hamiltonian formally written as Eq. ([20]) 
with star (or vertex) operators A s , each comprising the 6 
nearest-neighbors to a site, and planar plaquette opera- 
tors B p , each involving 4 spins (see Fig. [2]). This Hamil- 
tonian is basically a D = 3 Ising gauge theory (— J2 P B p ) 
augmented by the sum of all local symmetry generators 
(— ^2 S A s ). The constraints of Eq. ([23"]) get replaced by 



y[a s =i, n b p =i, 

s peCube 



(47) 



where Cube includes the six plaquettes which form the 
cube. Since no constraint couples the vertex and plaque- 
tte operators, the partition functions is simply 



Z 3 D — Z3D | s i n g g aU g e X Z\Y) | s i n g. 



(48) 



Let us first show that KE3Z? does show a thermo- 
dynamic phase transition at T > 0. From Eq. (|48p 
the free energy is given by F SD (j3) = F 3D | sing galIge (P) + 
Fid isin g (/3), where F 3D | sing galIge ((3) = F 3D \ S m g (f3*) with 
the dual (3* satisfying sinh2/3sinh2/3* = 1. Clearly, the 
free energy of KE3D displays a T = singularity coming 
from Fid ising (/?) and a finite-T singularity resulting from 
F 3 D Ising gauge(/3) at C = 0.761423. 

We now show that the D = 3 Ising Gauge theory dis- 
plays (at least) rank n = 8 TQO. [See the definition of 
rank-n TQO given in the Introduction.] Let us start by 
writing 8 GSs of the D = 3 Ising gauge theory on a cubic 
lattice of size L 3 which is endowed with periodic bound- 








1 





Cube 

FIG. 2: D = 3 extension of Kitaev's Toric code Model. On 
each bond (ij), there is an S — 1/2 degree of freedom indi- 
cated by a Pauli matrix Sij. Notice that elementary plaque- 
ttes are planar while stars involve 6 neighbors in the x, y, and 
z space directions. A Cube involves 6 plaquettes. 



ary conditions 



ce{+++} 
= ^ E l c >' 



= aa £ | 

ce{ — } 



(49) 



In Eqs. (|4"9"]) . the states {|c)} span all states in 
the afj basis which (i) have B p = 1 for each plaque- 
tte p [see the definition of B p in Eq. Q2ip] and (ii) 
lie in a specific topological sector and J\f is a uniform 
normalization constant. The eight topological sectors 

[(+ + +), (+ + -),•■■ , ( )] are labeled by three Toric 

invariants (fJ. — 1,2,3) 



n 



(50) 



Each of the cr^ eigenstates |c) is an eigenstate of 

the three Toric operators In Eq. (|50|) . each of 

the three cycles (+1,2.3 is a path which circumscribes one 
Toric cycle (e.g. a cycle along each of the three cubic 
axes). Each of the states in Eq. (|4"9"|) transforms as a 
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singlet under all of the cubic lattice star operations 

a s = n ( 5i ) 

z£star(.s) 

with the product above performed over all 6 bonds (si) 
which have the vertex s as one of their endpoints. This is 
so as the star operations of Eq. (|5Tj) link states |c) within 
the same topological sector. That is, if 

A s \c) = \c') (52) 

then, as [Z M , A s ] = for all fi and s, 

Z M \c) = Z M \J). (53) 

Given the invariance of the 8 GSs of Eq. (|4"9")) [and thus 
of any superpositions thereof] under all of the d = 
symmetries {^4 S }, we may proceed to demonstrate TQO. 
To this end, let us decompose any quasi-local operator 
V into the component invariant under all of the d = 
GLSs of Eq. (fSTj) (labeled by Va-o) and component Va-i. 
which does not transform as a singlet under all of these 
operations: V = Va-o + Va-±- For any state \g) which lies 
in the 8-dimensional space spanned by the states of Eq. 
([4"9"| . and which transforms as a singlet under all d = 
GLSs G, we have 

(g\VA;±\g) = 0. (54) 

All that we need to consider are thus the local (d = 0) 
symmetry invariant components Va-.o of the quasi-local 
operator V. The quasi-local operators invariant under all 
of the symmetries of Eq. (|5Tj) are built out of product of a 
finite number of operators {A s } and plaquette operators 
{Bp} (see Eq. (|2T|) '). We must now show that any such 
quasi-local operator Va-o attains the same expectation 
value in each of the states of Eq. (|4"9"|) . To this end, we 
consider the following three connecting operators 

X»= J] ofi- ( 55 ) 

These operators are the D = 3 extension of the operators 
X h2 of the D = 2 Z 2 gauge theory [Eq. pi])], 

In Eq. (|55|) . the product is taken over all bonds (ij) 
which lie in planes perpendicular to the cubic direction 
e M . These operators commute with one another and are 
symmetries of the Ising gauge Hamiltonian: [H, X M ] = 0. 
Acting with a particular on any state (e.g. any vortex- 
less a z eigenbasis state |c)) which is an eigenstate of all 
{Bp} operators with unit eigenvalue leads to states which 
are eigenvectors of {B p } with unit eigenvalue. More gen- 
erally, for any operators A s and B p , we have 

[A s ,X^] = [B p ,X )i ]=0. (56) 

As a consequence of Eq. (|55|) and the fact that all quasi- 
local operators Va-o are multinomials in {A s ,B p }, we 
have that 

[VA;0,X fl ]=Q. (57) 



Moreover, these operators satisfy the following algebra 
with respect to the Toric symmetries {Z^}: 

{A>, Z^} = , [A>, Z v ) = , n ^ v. (58) 

As a consequence of Eq. (|58|) , we see that the 8 states of 
Eq. (|4*9")) are related to one another by these operators. 
For instance, 

\g++-) = ^3b+++), 

!<?-+-) = X 1 X 3 \g +++ ), (59) 

etc.. Therefore, the 8 d = 1 GLS operators FJ M 
with rifj, = 0, 1 form a d — 2 group Q (of a Z2 x Z2 x 
Z2 character). These operators suffice to link all of the 
states of Eq. (|49|) with one another. By unitarity, these 
generators also link any set of 8 orthogonal states in the 
space spanned by Eq. (|49|) . 

By Eqs. ([541 IF7|) . the expectation value of any quasi- 
local operator V is the same in all GSs spanned by the 
n = 8 GSs of Eq. (H. Thus, the D = 3 Ising gauge 
theory exhibits (at least) rank-n, = 8 TQO. 

We note, in passing, that considerations similar to 
those above may be enacted for general Ik gauge theories 
(k = 2, 3, 4, • • • ) on a hypercubic Z?-dimensional lattice. 
This theory is defined by 

H ^ = -\ ^2(UijU jk U k iUu + h.c.) (60) 
p 

where, on every link (ij), a parallel transporter 

U iS = e w », 

% = 2-imij/k (61) 

with nij an arbitrary integer. The elements Uij of Eq. 
(|6"T) satisfy a Z& algebra. Here, instead of Eq. (fB"5")) . we 

set 

n ei ^ Lti - ( 62 ) 

(u)J-e fl 

In Eq. (|62|) . is the generator of rotation of the Z& 
variable on bond (ij). Replicating the proof given above, 
we find that the Z& gauge theory on the Z?-dimcnsional 
lattice exhibits (at least) rank-n = k D TQO. 

It remains to prove that at any finite T > no topo- 
logical symmetry operator may acquire a non-vanishing 
expectation value. For example, 

<^>=< n 4)=°' ( 63 ) 

for {Cfj,} loops around the Toric cycles. [Equation (|63l) is 
a particular realization of Eq. (JT2J) .] Due to the decou- 
pling of the plaquette ({B p }) and vertex ({A s }) opera- 
tors, the expectation value (Zc 1 Zc 2 ) is given by its value 
for a classical D = 3 Ising gauge theory. However, as seen 
by large and small coupling expansions [l6j . the correla- 
tor (Z Cl Z C2 ) = (U.{ i j)eC 1 a ijU(ij)eC 2 a ij) vanishes in 
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the D = 3 Ising gauge theory as the bounding contours 
Ci t 2 are taken to be infinite. [For finite size systems, 
(ZcJ) vanishes as no SSB is possible.] In (i) the confined 
phase [(3 < (3 C 2\ this correlator vanishes exponentially in 
the area of the minimal surface R bounded by the Toric 
cycles Ci,2 while in (ii) the deconfined phase \j3 > /3 C 2], 
the pair correlator (Zp ± Zp 2 ) vanishes exponentially in 
the total length of the contours C\ and C^- In the limit 
of far separated contours C\ and C2, both (i) and (ii) 
reaffirm Eq. (|6"3")) . Equation (f6"3"| is suggestive of a finite 
autocorrelation time at all positive temperatures. 



C. Kitaev's Honeycomb Model 

Kitaev's model on the honeycomb lattice [f| is defined 
by the following S = 1/2 Hamiltonian (Fig. [3]) 

H Kh =- J x ]T Of Of - Jy -J*J2 ( 64 ) 

a: — bonds y— bonds z — bonds 
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FIG. 3: Kitaev's model on a honeycomb lattice and three 
types of bonds. On each vertex there is an S = 1/2 degree 
of freedom indicated by a Pauli matrix Sj. C represents an 
arbitrary contour drawn on the lattice. 

Here, we find that Elitzur's theorem mandates that all 
non- vanishing expectation values must be of the form 25 1 



Or 



(65) 



with C any contour (set of contours) drawn on the lattice 
and 7 is the bond direction which is orthogonal to the 
path C. Let us consider embedding this system on a 
torus with g handles (a torus of genus g). For closed 
contours C which do not span an entire Toric cycle, Oc 
is not independent of the code's stabilizer. In fact, there 



is a rather simple connection between these operators 
and the local anyon charge which measures the number 
of local high energy defects. For a closed contour C which 
is an elementary closed hexagonal loop, we have Oh which 
is the anyon charge associated with a given hexagon h. 
Rather specifically, 



O h = aloZolalofol 



(66) 



with 1 — 6 labeling the vertices of any given hexagon of 
Fig. El The spin polarization directions fi for all spins 
in the product of Eq. (|6"6")) have been chosen to corre- 
spond to the single bond direction (i (x, y, or z) that is 
attached to the site i and does not lie on the hexagonal 
path. 

As shown by Kitaev, within the GS sector, (Oh) = 1 for 
all h. Extending the arguments of the strong and weak 
coupling expansions of gauge theories [l6| . we now find 
that the correlator between any two contours (OcOc) 
scales as (i) e~ ClA at high temperatures with A the area 
bounded by C and C and c\ a positive constant and, if 
an ordered phase exists, scales at low temperatures as (ii) 
e -c 2 (\c\+\c'\) witn |C| + |C"| the total perimeter of the two 
loops and C2 another positive constant. Thus, for cycles 
C which span the entire lattice we find as before that no 
d = 1 loops of the form of Eq. (frJS"]) can attain a finite 
expectation value. 

Let us now consider closed contours C1.2 which span 
the same Toric cycle. For a system with periodic bound- 
ary conditions, at T = 0, we now have for closed contours 

Cl,2, 



<n°? n °i> = <n d *> = l 

ieCi jec 2 heD 



(67) 



Here, D is the domain bounded by two contours 
which go around an entire cycle of the torus [see e.g. the 
thick solid line in Fig. [3], Taking the separation between 
the two contours C\ and C2 to be very large, we have 
that, at T = 0, 



<n°? n "?>->< n°f>< ih: 

ieCi jec 2 ieCi jec 2 

This implies that, at T = 0, 



(68) 



(69) 



At T > 0, no SSB may occur (also in the thermodynamic 
limit, see Eq. (fTg)) ) and all quantities of the form Eq. (|6"9")l 
must be zero. 

We now write down a formal finite temperature so- 
lution to the partition function and comment on low- 
and high-T series expansions. Using fermionization [25[, 
Kitaev's Honeycomb model can be cast as a model for 
fermions {d r } on a square lattice with a site-dependent 
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chemical potential i] r = ±1. 

H Kh = J X J2 ( d r + d r) { d t+e x ~ d r+e x ) 
r 

+ JyJ2( d l+ d r) ( d t+e v - d r+e y ) 
r 

+ J z J2Vr(Hd r -l) 

r 

= <tfM n d + dtW n dt + dW}d. (70) 

The fermionization of pBj has been recently invoked to 
attain very interesting results in extensions of Kitaev's 
Honeycomb model [2g. In Eq. ([TO")) , {r} denote the 
centers of the vertical bonds. The unit vector e y connects 
two z-bonds and crosses a y-bond, see Fig. [3] A similar 
definition holds for e x [25[ . The form of Eq. (f70|) is very 
similar to the Fermi representation of the D = 2 Ising 
model. In the last line of Eq. l[70"]). 

dt = ({4», d=({d r }). (71) 

With A = (dt,d) [13 and 

N _( M v 2W n \ 

the partition function for the Fermi bilinear of Eq. (1701) 
is easy to write down 

Z Kh = J2 e" ftr [M " ] | det(e^" + 1)| 1/2 . (73) 

{>?,•} 

Within the GS sector (the one with no anyons), the sys- 
tem is translationally invariant and its spectrum [f| is 
that of a p-wave type BCS pairing problem [25L l26j|. 

Due to the sum over all of the 2 N "' 2 Ising configura- 
tions {r) r }, the complete partition function of Eq. (|73|) 
is non-trivial. The existence of a transition as temper- 
ature is varied is not as immediate as in Kitaev's Toric 
code model. Translational invariance appears only for 
uniform r\ r = 1 or — 1 for all r and by a simple unitary 
transformation to a system in which r) r is constant on 
entire horizontal lines (25[. Nevertheless, bounds on cor- 
relators are easily established: any fermionic correlator 
computed with the full partition function of Eq. ((73"]) is 
bounded from above by its value when computed within 
the sector {r] r } which maximizes its value. The finite T 
fixed sector correlators for the quadratic Hamiltonian of 
Eq. ([70)1 can be computed with the aid of Wick's the- 
orem. If the GS sector (that with rj r = 1 everywhere) 
is gapped then at all temperatures, the correlation func- 
tions exhibit exponentially decaying correlations. 

We may expand Eq. (|73[) in low- and high-T se- 
ries. A high-T series may be derived by an expansion of 
| det(e /3Ar " + 1)| in powers of (3. At low temperatures, the 
GS terms correspond to the phase T] r — 1 for all r which 
just reproduces the p-wave type BCS result of 0, [25|. 
At finite temperatures, we allow for r\ r ^ 1. The lowest 



energy terms correspond to a few vortex pairs [r\ r = —1) 
which are tightly bound. We may diagonalize the Hamil- 
tonian to determine the spectrum for these vortex pat- 
terns and find the low-T corrections to the GS results. 
In such a manner, we may write down a low-T series 
expansion for the correlation function between fcrmion 
pairs (which corresponds to string correlators in terms of 
the spin variables of Kitaev's model) . The explicit high- 
and low-T series expansions will be provided elsewhere. 
For now, let us note that within the gapped phase, the 
system exhibits a finite correlation length for these string 
correlators [2{| . It is clear that if a low-T series expansion 
about the GSs is possible, then at sufficiently low tem- 
peratures, corrections to the T = result for the string 
correlators can be made arbitrarily small. In particular, 
when such an expansion is carried out about the GSs of 
the T = gapped phase, the string type correlators of 
[l§] still display a finite correlation length. [As we dis- 
cussed above, that these string correlators must indeed 
display a finite correlation length also follows from the 
form of the fermionic correlator in the {r] r } sector which 
maximizes its value; this value of the correlator provides 
an upper bound on the correlation function of the com- 
plete (unpinned {r) r }) system]. The correlation function 
on the righthand side of Eq. is related to a lim- 

iting form of such string correlators in which the string 
(defined here by the contour C\ ) wraps around an entire 
Toric cycle. 

In what follows, we comment on the prospect of finite- 
T phase transitions in Kitaev's Honeycomb model as it 
follows from our expansions. As we show in the Ap- 
pendix, in the absence of a phase transition, the topolog- 
ical quantum error detection condition of Eq. (JSJ) [and 
its equivalent Eq. (|A1[) ]. cannot be satisfied. If no tran- 
sitions exist, then hardware fault-tolerance at finite tem- 
peratures may be impossible to achieve. The T = gap- 
less phase - the phase which is not robust - seems to 
be amenable to a phase transition as the temperature of 
the system is increased. In the gapless phase, a transi- 
tion from a system with oscillatory algebraic correlation 
between fermions at T = to one with exponential cor- 
relations at high temperatures is suggested by the low- 
and high-T expansions. For states deep inside the gap- 
less region [H, [25| a low-T expansion (that includes sec- 
tors with a non-vanishing number of anyons) suggests 
that low-energy states may still carry a vanishing gap 
and lead to a T > phase with algebraic correlations. 
Thus, a finite-T transition between a high-T phase with 
exponentially decaying correlations between fermions to 
a low-T stripe type phase with oscillatory algebraic cor- 
relations may exist. 

In the absence of applied external fields which intro- 
duce a gap, the gapless sector of the theory is, however, 
unstable to perturbations. Of greater promise - insofar 
as stability to local perturbations - is the gapfull region of 
the Kitaev model. Regretfully, no transition of the type 
found for the gapless phase is evident for the gapfull sec- 
tor of the Honeycomb model. Thus, albeit satisfying the 
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T = condition of Eq. (f4]), in the absence of a finite- 
T phase transition, the condition of Eq. ([5]) cannot be 
satisfied at any finite temperature. 

Our results on the low- and high-T expansions, which 
suggest that the gapfull sector of the Kitaev model does 
not exhibit a finite-T phase transition, can be further for- 
tified and linked to expansions carried by earlier works. 
These expansions elucidate how the gapfull phase will 
not break ergodicity and cannot achieve hardware fault- 
tolerance. The correspondence to our earlier results con- 
cerning Kitaev's Toric code model is made clear by a sim- 
ple mapping. In the parameter region J z J x ,y, & per- 
turbative expansion gives the following effective Hamil- 
tonian || 



H, 



Jll 



eff 



16J? 



(74) 



In other words, we recover Kitaev's Toric code model 
of Eq. (|20|) . As a result, we have as before that the 
logical operators (Xi 2 and Z\^) have a finite (system- 
size independent) autocorrelation time r [Eq. (139|) ]. 



IV. THERMALLY STABLE ORDERS IN THE 
CASE OF HIGH-DIMENSIONAL BROKEN 
SYMMETRIES - A CASE OF SIMPLE NON-TQO 

Some systems with non-local symmetries exhibit a 
robustness to thermal fluctuations. These systems are 
natural contenders for satisfying the conditions of Eq. 
©. In all of the cases that we examined, however, ther- 
mal fragility still appears. 

Let us start with the simple D = 2 orbital compass 
model. Its Hamiltonian is given by 



H 



3 



3 " j+e* 



(75) 



which emulates the direction dependent interactions re- 
sulting from the anisotropy of electronic orbitals. In this 
model, exchange interactions involving the x component 
of the spin occur only along the spatial x-dircction of 
the lattice. Similar spatial direction dependent spin ex- 
change interactions appear for the z components of the 
spin. Apart from a global reflection which only appears 
for the isotropic point J x = J z (and which may be bro- 
ken at low T), the anisotropic orbital compass model 
has the following symmetry operators O m = JljeCj, 
for n = x, z. denotes any line orthogonal to the e M 
axis. As O m involves 0(L X ) sites, they constitute d = 1 
symmetries. The spectrum of the D = 2 orbital com- 
pass model is gapless [28[ in the thermodynamic limit 
due to the existence of the symmetries that it possesses 
@j EH [29[ . What is important for relaxation processes 
(much as was evident in the analysis of Kitaev's model) 
is not at all the energy spectrum but rather the free en- 
ergy. It was found that the large S renditions of the 



isotropic D = 2, 3-dimcnsional orbital compass model, 
although highly degenerate in the low-energy spectra in 
the thermodynamic limit, exhibits free energy barriers. 
An entropic stabilization sets in which at sufficiently low 
yet finite < T < T* leads to a sharp phase transition 
with a divergent relaxation time (30j . It is noteworthy, 
however, that at the isotropic point, the system is not 
topologically ordered as a global inversion symmetry 



O 



Reflection 



(76) 



can be broken. It is this very symmetry which is bro- 
ken below T* . Here, the local nematic order parameter 
operator 



3 (J 3+e I 



(77) 



attains a finite expectation value at low T |6j, [l_2j, |29|, [3 



(N 3 ) a ? (N 3 ) Pl T < T* 



(78) 



Thus, measurements of this quasi-local order parameter 
within the GS sector lead to different answers for dif- 
ferent contending low-energy Gibbs states (labeled by a 
and /?). Thus, Nj is not a contending logical operator - 
it does not satisfy Eq. (jTJJJ) . For all non-local operators 
M we have that at any finite T, (0 M ) = on any lattice 
(whether it is finite and no SSB is possible or on the in- 
finite size lattice where the perimeter law scaling implies 
a vanishing = 0). 

Similar arguments apply for the D = 3 orbital com- 
pass model, meaning that the d = 2 symmetries (as well 
as a three-dimensional reflection symmetry) will render 
the TQO non-existent at finite temperatures, although 
the relaxation time is obviously divergent. We recently 
became aware of a paper by Bacon [31| where he pro- 
vides evidence (by using a mean-field approach) of a self- 
correcting system in a D = 3 relative of the orbital com- 
pass model. The reason why this system is claimed to 
be a high-temperature self-correcting quantum memory 
(i.e. one without additional quantum error correction) 
is because, at the mean-field level, low-energy excitation 
energies may scale with the perimeter of the defected do- 
main. As in the D = 3 orbital compass model, Bacon's 
model display d = 2 (discrete) g?-GLSs. The problem we 
see with this argument is that these symmetries can be 
broken giving rise to an ordered phase (i.e. one with a 
Landau order parameter). Then, condition Q is going 
to be violated and no TQO will be present. As in other 
cases in which symmetries can be broken, here Eq. (JTTJJ) 
- a version of which is related to the weaker version [Eq. 
©] of the finite-T error detection condition of Eq. ([5]) - 
will be violated. 



V. CAN THERMAL FRAGILITY BE 
DEFEATED? 

It is clear that a physical system that is in an ordered 
thermodynamic phase, i.e. having a non-vanishing Lan- 
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dau order parameter (i.e. breaking a global symmetry), 
displays some sort of robustness against local errors due 
to the collective nature of the resulting order. This prop- 
erty is used to build robust classical memories. We have 
seen that topological quantum memories seem to be ther- 
mally fragile, therefore, it seems plausible that by build- 
ing systems that display an order parameter we could 
defeat it and obtain robust quantum memories. In gen- 
eral, there seems to be a catch regarding stable quantum 
topological memories. According to Eq. (JSJ), in order 
to ensure error detection of local perturbations at finite 
temperatures, the system must display at least one fi- 
nite temperature phase transition (see Appendix |A"|) . We 
generally find that (2),) =0. As a matter of princi- 
ple, even if we allow for operators {R a } of a topological 
character (that is, non-local operators) which are not in- 
dependent of the code's stabilizer, unless a low-T series 
expansion about the T — state is void, we will always 
get a perimeter-law type scaling for the quantum topo- 
logical observables, see e.g. Eq. ([T8|). This expansion 
for local Hamiltonians [see e.g. Ref. [16[] always gives us 
the result that the topological observables vanish for large 
loops (for which local errors can be avoided at T = 0). 
The expansion about T = can be void if there is a T = 
transition on an infinite size system. Such an occurrence 
matters worse once again as it implies that the robust- 
ness at T = does not obviously imply a robustness at 
any finite temperature (even at temperatures T = + ). 
On any finite size system, there is no SSB and as before 
(Tft) = identically. We illustrated that the same also 
occurs in the thermodynamic limit of these systems. Er- 
godicity cannot be broken insofar as the logical operators 
{T^} can detect it. 

We do not see an obvious way to avoid this conun- 
drum both for the viable logical operators {T^} which are 
independent of the code's stabilizer and more generally 
also for general non-local operators {R a } that can encode 
topological operations. It should be emphasized that all 
of our results pertain to rigorous bounds. Although finite 
(and size-independent) autocorrelation times seem to be 
mandated, if situations exist in which A 3> fceT may 
be realized then the autocorrelation times can be made 
large (see e.g. the explicit form of Eq. P0"| ). Based on 
heuristics and static arguments, a recent review provides 
estimates on these ratios in some candidate systems [32j . 

VI. SUMMARY 

In the present article we expanded on the concept of 
thermal fragility and applied it to several case examples 
including Kitacv's Honeycomb model and extensions to 
higher spatial dimensions. Our main results are 

(i) We introduced and discussed quantum error detec- 
tion at finite temperature as it applies to topologically 
ordered systems [Eq. ©]. 

(ii) We find that a system cannot satisfy the finite tem- 
perature error detection criteria of Eq. (jSJ) unless it dis- 



plays a finite temperature phase transition. It is impor- 
tant to emphasize that the existence of a finite tempera- 
ture transition (or several finite temperature transitions) 
is a necessary but not sufficient condition for Eq. (JSJ) to 
hold. 

(iii) In general TQO systems, whether they display 
finite-T phase transitions or not, and whether they have 
a finite size or are systems taken in the thermodynamic 
limit, all non-trivial topological operators may have a 
vanishing expectation value at all non-zero temperatures: 
{Tfj} = 0. This suggests that the code might self-correct 
only up to finite cutoff times that are governed by tem- 
perature dependent effects. Depending on the relative 
size of the gaps to the temperature, this autocorrela- 
tion time can practically be made very large. This 
may go hand in hand with (system size independent) 
finite bounds on the autocorrelation times r as adduced 
from the values of 1(7^ (0)2^(^)1. For times \t\ > t, 
1(^(0)2^(^)1 tends to zero [Eqs. ([151 1551 1301 HS))] . In- 
formation stored in the correlators of {!),} may be lost. 
Insofar as the non-local topological operators {7^} can 
detect, the system becomes ergodic at large time. Al- 
though of, apparently, less pertinence to ideal quantum 
memories, the same result is found for general operators 
{R a } which encode (non-local) topological operations on 
infinite-size systems but are not necessarily independent 
of the code's stabilizer, Eq. (p~8j) . or correspond to sym- 
metries. 

(iv) Kitaev's Toric code model can be solved exactly. 
At all non-zero temperatures, the autocorrelation func- 
tion exhibits a finite equilibration beyond which ergodic- 
ity sets in and all information is lost. The explicit forms 
for the autocorrelation functions of Kitaev's Toric code 
model are given in Eqs. (|39l l40l |43| . 

(v) Kitaev's Honeycomb model can formally be solved 
at finite temperatures. The main utility of this formal so- 
lution is that it enables low- and high-temperature series 
expansions. These expansions suggest that this model 
may not be stable to thermal fluctuations (in the sense 
ofEq. ©). 

(vi) We showed that some candidate systems which are 
stable to thermal fluctuations often display local orders. 
These local orders render the system unstable to errors 
from local perturbations and the topological protection is 
lost. 

A vexing question is how to quantify topological quan- 
tum error conditions at finite temperatures. The con- 
nection between rigorous bounds on quantum error de- 
tection and the more physical (although rigorous) results 
presented in our work is far from clear. In our opin- 
ion, this question needs to be quantitatively addressed. 
Given the physical results that we derived in our work, 
it is unclear to us if there is no way to overcome thermal 
fluctuations and still maintain quantum self-correction 
that is immune to local perturbations. We hope that 
our work will stimulate further studies which will ad- 
dress many related questions: Can d — 2 GLSs lead to 
robust quantum memories? What is the relation between 
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a finite lifetime r and the vanishing expectation value of 
the logical operators T M ? What happens when a certain 
amount of disorder is present and thus, in general, d- 
GLSs are not (exactly) present? In that case it is not 
true that for a finite system the expectation value of the 
logical operators vanishes since, for example for Kitaev's 
Toric code model [Hk + V, Z M ] ^ where V represents 
the disorder. We conjecture that, nevertheless, in the 
thermodynamic limit those expectation values vanish, for 
example (Zu) — > 0, and more importantly the autocor- 
relation times can generally remain finite as long as V 
is an arbitrary quasi-local operator (and the amount of 
disorder is not big). Indeed, in our earlier work [||, we 
showed that system's thermodynamic behavior and re- 
sponse functions may remain adiabatic even when their 
rf-GLSs may be broken. 

The main message is that we need to be careful when 
designing a topological quantum memory. We are forced 
to consider finite temperature effects since a real system 
is subjected to temperature effects. The common lore 
that a finite gap can indefinitely protect states (and thus 
quantum information) is not always correct. The limit- 
ing size of the autocorrelation time can often be system 
independent. If very low t emp eratures may indeed be 
achieved relative to the gap [32| then the autocorrelation 
time cutoff (e.g. Eq. (j4"0)) ') can be made large (albeit still 
finite). It may also be that systems with a well crafted 
environment can be constructed; such environments will 
disallow rapid thermal equilibration and thus keep the 
autocorrelation times high. Possible realizations of such 
environments can include a frequency mismatch between 
the natural resonant frequencies of the TOO system and 
the environment in which it is embedded [331 ] . 

In some of the Appendices that follow, we elaborate 
on points related to the main results reported here. In 
Appendix [XJ we prove that in the absence of a phase 
transition (a singularity in the free energy), the finite 
temperature error detection condition of Eq. (|8|) might 
not hold. In Appcndix[Bj we discuss a trivial consequence 
of the Licb- Robinson bounds to topological quantities. In 
Appendix [Cj we sketch how rigorous bounds for the ex- 
citation of anyons may be derived in some disorder free 
systems in thermal equilibrium. The derivation outlined 
here is intimately linked to Kitaev's Honeycomb model of 
Section HTlCl However, the result is far more general. In 
Appendix [Dl we discuss thermal equilibration time (the 
time for the system to become ergodic) in cases when the 
system displays or is close to a critical point and general 
scaling laws become apparent. 

Note added in proof. After completing the current 
work and much after [6( , we became aware of the inter- 
esting independent work of Ref. [2(| which raised con- 
cerns about the reliable storage of quantum information 
at finite temperature. Reference [20| invoked the KMS 
equations in the study of Kitaev's Toric code model to 
independently also arrive at and fortify one of the con- 
clusions of Ref. Q (that of the last line of Eq. (| 3 1 [) ) from 
a different approach. 
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APPENDIX A: TOPOLOGICAL QUANTUM 
ERROR DETECTION AT FINITE 
TEMPERATURES 

We show that any system that does not exhibit a fi- 
nite temperature (T > 0) phase transition (regardless 
of whether it exhibits TQO or not) necessarily violates 
the quantum error detection condition of Eq. ([5]) at all 
temperatures T > 0. In other words, for all T > 

W=[p x ' 2 Vp l / 2 ,%]^Q. (Al) 

Here, we recall that p = Z^ 1 exp[— H/{k,BT)\ with Z = 
tr [cxp[— (3H]\ = cxp[— [3F] with F the free energy. In 
the absence of a phase transition, the free energy may be 
expanded about its high temperature (small /?) limit 

oo 

F = Y J F k p k . (A2) 

fc=0 

All terms in the series of Eq. (|A2[) are c-numbers. If the 
system displays a phase transition the expansion of Eq. 
(|A2|) has a finite radius of convergence ((3 C ). If the system 
displays no phase transition, the series of Eq. (|A2[) is 
everywhere convergent. 

The proof of the assertion of Eq. (|A1|) for the case a 
divergent [3 C (a system with no finite temperature tran- 
sitions) is trivial. For any finite system at temperatures 
T > 0, the density matrix p can be expanded to all orders 
in /3 = l/{k B T). We have that (p = exp[-/3(H - F)}) 

W= [V, %] - !/?[{(# - F ), V}, f M ] + ^[{F U V}, %] 

+ ^-{{(H- F ),{(H F ), U}}, T„] + 0(P 3 ). (A3) 

At any fixed chosen (3 for which F(/3) is non-singular, the 
expansion of exp[— (3(H — F)] in powers of (3, with F{(3) 
held at its fixed value at the chosen (3, has an infinite 
radius of convergence. In Eq. (|A3j) . we collect terms of a 
fixed power of (3 when we use the expansion of Eq. (|A2[) . 
In the absence of a finite temperature phase transition, 
the series of Eq. (|A3[) converges for all (3. If W van- 
ishes on a dense set of points along the real axis with 
arbitrary then as W is an analytic function then it must 
also vanish everywhere in the complex /3-plane. All of 
the derivatives of W and, in particular, its zeroth order 
value Wo = [V, T M ] must vanish. However, Wq is inde- 
pendent of the system in question (it is independent of 
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H, although [H, T) t ] = 0). For any non-trivial algebra 
captured by the operators {T^}, there exists quasi-local 
operators V that (i) have their support on the same re- 
gion where {T^} are defined, and (ii) do not commute 
with {T^}. Eq. (|A1[) docs not preclude a commutator 
which vanishes in the thermodynamic limit. The im- 
possibility of Eq. ijSJ) in the thermodynamic limit in all 
systems which exhibit SSB detectable by quasi-local op- 
erators follows from Eq. @. In all finite size systems, 
there is no SSB (and the free energy is analytic - f3 c = oo) 
and Eq. (|A1|) cannot be satisfied at any non-zero tem- 
perature. A similar conclusion holds for systems which, 
even in their thermodynamic limit, do not exhibit T > 
transitions. Kitaev's Toric code model is precisely such 
a system. It is important to emphasize that the appear- 
ance of one (or more) T > transitions is a necessary 
but not sufficient condition for Eq. §8§ to hold. 

APPENDIX B: EXTENSIONS OF THE 
LIEB-ROBINSON BOUND TO NON-LOCAL 

TOPOLOGICAL QUANTITIES 

Locality in quantum spin systems is often well cap- 
tured by a bound on the commutator of local obscrvables 
with disjoint supports. Early on, Lieb and Robinson pro- 
vided estimates and bounds on precisely such commuta- 
tors on nearest-neighbor spin systems on general graphs. 
These naturally provide bounds for the speed of propa- 
gation of quantum information. Additional works elab- 
orated on the initial findings of Lieb and Robinson and 
examined extensions to several bosonic and fermionic sys- 
tems QJ]. 

Stated formally, the Lieb- Robinson bound asserts that 
given two local operators V a and V b in different regions (a 
and b), the operator norm of the commutator in gapped 
systems 

l|[Va(o),H(t)]|| <A(yo\ min )\\v a \\\\v b \\e-^- v ^. (Bl) 

In Eq. (|B1[) , / is the distance of the shortest path linking 
regions a and b (more precisely the number of links on 
such a path), vol m i„ = min{|a|, \b\} denotes the number 
of vertices in the smaller of the two sets a and b. The con- 
stants A, v, and £ depend on the maximal coordination 
of a given vertex of the lattice and the largest local terms 
max{||/ijj||} which appear in the argument of a nearest 
neighbor Hamiltonian system 

# = (B2) 

Similar results are found to pertain to finite temperature 
correlators for systems with local Hamiltonians. 

The extension of Eq. (|B1[) to non-local quantities is 
straightforward. Let us consider the general product of 
local operators over an L d -dimcnsional volume 

f„=l[Va. (B3) 



If the operator norm of all Vs is bounded from above by 
unity then the operator norm of the commutator 

||[ri(o),r 2 (t)]|| < I^RiII-Ral , max \\[V a ,(0),V b ,(t)}\\ 

< AlExH^KvoU^llK'll || ||e- C'-^— (B4) 

with a' € i?i and b 1 £ R 2 . In Eq. (|B4[) . Z min denotes the 
smallest distance between two points - one which belongs 
in i?i and the other in R 2 . The first line of Eq. (jB4)) 
follows from a long-hand form of the expansion of the 
commutator. As T\ contains |i?i| terms and Tg has |i?2| 
terms, there are a total of |i?i||i?2| terms each containing 
a string of length (|J?i| + I-R2I — 2) of local operators 
tya\ which bracket a single commutator between two 
local operators in the disjoint regions [V^, Vb']. If the 
operator norm of all Vs is bounded from above by unity 
then the bound of the last line of Eq. (|B4[) follows from 
Eq. flm]) . 

As Z m i n — > 00, the equal time commutator between two 
distant topological operators vanishes. In some cases, the 
speed of propagation v could be related to the velocity 
of defect motion (e.g. domain walls). 

We can similarly bound the operator norm of the com- 
mutator which appears in Eq. (jSJ). Here, we simply have 
that 

\\[f»(0),V b (t)]\\ < \R\ max \\[V a ,(0),V b (t)}\\ 

< A\R\(vo\ min )\\V a <\\ \\V b \\ e -«— (B5) 

Although the influence of events far away from the sup- 
port of T M (region i?i) can become nil (as seen from Eq. 
(|B5|) . the cumulative effect of fluctuations on and prox- 
imate to Ri can (and indeed generally does) lead to a 
vanishing (X^). 

APPENDIX C: THE PROBABILITY OF 
OBSERVING ANYONS 

We now sketch how rigorous bounds on the probability 
of observing anyons may be derived in one of the models 
that we discussed in this work (Section IIII CP - Kitaev's 
Honeycomb model Q . Although some of the expressions 
given below are special, the result seems to be of greater 
applicability. 

As shown by Kitaev, the GS configuration is the anyon- 
free configuration in which (Oh) = 1 for all plaquettes 
(hexagons) h. Oh is the product of the two rj r quanti- 
ties which live on the center of two vertical bonds which 
appear in h |25l[. As (i) this system exhibits reflection 
positivity (RP) [3J| for the Ising-type free energy (-Fising) 
in the argument of Eq. (|73l) , 

e -PF ls ; ne [ M ] = e -|tr det ( e 0AT„ + !)|l/2 (C1) 

(ii) the minimum of the spectrum of an L x L slab of the 
system [of the square lattice of [2^|] with open boundary 
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conditions within a given topological sector {Oh} is dis- 
crete, and as (iii) within the GS the system is anyon free, 
there must exist a length L c such that in all slabs of size 
L > L c the lowest free energy attainable over all states 
which are not anyon free is larger by a finite gap Al > 
than the free energy within the anyon free infinite size 
system endowed with periodic boundary conditions. As, 
for example, in Ref. (3q and in particular [36| . this im- 
plies that the probability of finding anyons (or vortices in 
36] ) is exponentially suppressed at sufficiently low tem- 
peratures. This is done by tiling the lattice with all L x L 
blocks and noting that the probability for a given config- 
uration is bounded from above by 



p = exp 



N r . 



(C2) 



where w is a constant of order unity, iV& is the number 
of bad blocks which contain (at least) one anyon and N c 
is the number of L x L blocks to which a given bond 
is common. Combined with RP of the Ising-type Gibbs 
weight of Eq. (|C1[) . we can prove that the probability 
of anyons is exponentially damped in f3. Thus, although 
a gap is generally required for quantum protection from 
perturbations it may render the probability of excitations 
of such anyons exponentially small. Equation (|C2[) is a 
rigorous result for any uniform RP system [including that 
of Eq. (|64p]. It is noteworthy that the introduction of 
disorder may introduce anyons (and general defects) even 
at low temperatures. This has implications for impurities 
in Quantum Hall systems. 



APPENDIX D: THERMAL EQUILIBRATION 
TIME IN SYSTEMS WITH DISCRETE 
SYMMETRIES 

We can use the following heuristic metastability argu- 
ment to address the important question: How long does 
it take for a TQO system with d-GLSs which has been 
prepared in the GS (protected subspace) to equilibrate at 
a temperature T in the presence of d-dimcnsional defects? 
Stated alternatively, we now address the following ques- 
tion: how long can the system retain information before 
it is lost due to ergodicity? The general answer to this 
question is quite complex. One can use qualitative ideas 
borrowed from the theory of mctastable states, where the 
thermal equilibrium state is the global minimum of the 
free energy F while the GS manifold is the suitable sub- 
space which does not minimize F but is supposed to be 
locally stable. The lifetime of the protected subspace r is 
clearly determined by the dynamics of the system when 
coupled to a thermal reservoir. The basic problem then 
reduces to characterizing the free energy barrier AF that 
obstructs the escape from the GS manifold. The crucial 
observation is that for a system of size L x L ■ ■ ■ x L for 
a discrete symmetry of dimension d > 1, the free energy 
barrier follows the form 



with a a positive constant and a the tension. The point 
is that the presence of low-dimensional symmetries lower 
the height of the barrier because of dimensional reduc- 
tion. The case d = 1 was considered in Section IIIII (al- 
though the energy cost was finite in that case the entropy 
scaled logarithmically with the system size). 

Let us assume that the probability to occupy an eigen- 
state \(f>j) of H } 



P 7 (t) = falUtpUfa 



(D2) 



with p the initial state and U a unitary evolution that in- 
cludes the coupling to the thermal bath, follows a master 
equation 



dP 7 (t) 
dt 



: ^(H/ 77 ,P 7 ,(t) 



W 7 / 7 P 7 (i)), (D3) 



with transition probabilities W 77 / satisfying the detailed 
balance condition (with equilibrium in a canonical en- 
semble) 



W 11 > exp[— {3Ey] = W-y'-y exp[— /3-E 7 ]. 



(D4) 



It is clear that the whole complexity of the dynamical 
problem is in the functional form of the transition prob- 
abilities. Let us assume that there is only one possible 
channel for the system to escape from the local minima. 

Kitaev's model is quite special in that it displays a 
T = critical point. What can we say about systems 
that have a finite transition temperature? 

Typically, in systems with a finite critical temperature, 
the equilibration time follows the Arrhcnius form 



exp[(3AF] 



(D5) 



with AF the free energy barrier for the defects. This, 
along with Eq. (|Dip leads to the usual form 



(D6) 



t - cxp[a/3L d -V] 



In what follows, we address forms of this equilibration 
time and the surface tension which governs it in several 
regimes. 



1. Equilibration time scaling near a critical point 

The considerations presented below apply general scal- 
ing arguments (see e.g. [13]) to general systems with dis- 
crete symmetries. For such systems, next to the critical 
point, 



a = k B T c /£ 



d-l 



(D7) 



with £ the correlation length. Near the critical tempera- 
ture T c the equilibration time scales as 



AF = aL 



d-l. 



(Dl) 



cxp[c(L/0 



d-li 



(D8) 
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with c another positive constant. We can insert the scal- 
ing of £ with the reduced temperature t = (T — T c )/T c 
near T c 

Z~\t\- V (D9) 
to obtain the relaxation time near the critical point 

r~ Aexv[cL d - l \t\ v{d -\ (D10) 

a. Finite surface tension when T — > 

In most systems (not all) having a finite T c , the surface 
tension tends to a finite value as T — > 0. Consequently, 
Eqs. (|Dll ID5[) imply that the equilibration time r is 
strongly divergent in both the system size L (fixed tem- 
perature) and the inverse temperature (3 (at fixed L). 



b. Temperature dependent surface tension 

A temperature dependent surface tension may appear 
in the orbital compass model and other systems which 
exhibit an order-out-of-disorder phenomenon (30l . [351 ] . 
Here, we can prove that at sufficiently low temperatures 



f3a > q > 0. (Dll) 



This, then, enables the proof of a finite-T phase transi- 
tion by the usual Peierls argument. However, although 
the relaxation time r is divergent in system size it does 
not, for a finite fixed L, diverge as the temperature 
T ^ 0. 
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